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Compact formulas for trapped-particle and passing-particle guiding-center orbits in axisymmetric 
tokamak geometry are given in terms of the Jacobi elliptic functions and complete elliptic integrals. 
These formulas can find applications in bounce-center kinetic theory as well as neoclassical transport 
theory. 



I. INTRODUCTION 

The presence of trapped-particle and passing-particle 
guiding-center orbits in axisymmetric tokamak geometry 
has had a significant impact on magnetic fusion research 
(see Refs. [20 for a brief historical survey) . Because of 
their fundamental importance in understanding neoclas- 
sical and anomalous transport in tokamak plasmas [3] , it 
is desirable to find compact analytic representations for 
the trapped-particle and passing-particle guiding-center 
orbits in axisymmetric tokamak geometry. In particu- 
lar, compact expressions might allow explicit calculations 
of bounce-angle averages in bounce-center kinetic theory 
beyond the deeply-trapped approximation. In ad- 
dition, a bounce-center Fokker-Planck collision operator 
might be constructed from the guiding-center Fokker- 
Planck collision operator Q, which would enable a ki- 
netic description of neoclassical transport beyond the 
zero-banana- width limit Q . 

Previous analytic representations have used a simpli- 
fied tokamak geometry [9| in which the magnetic sur- 
faces are unshifted concentric circular toroidal surfaces 
centered on the magnetic axis (located at a major radius 
R for the axis of symmetry). In this simplified tokamak 
geometry, trapped-particle and passing-particle guiding- 
center orbits [y] are represented in terms of complete and 
incomplete elliptic integrals [10| . The purpose of this 
work is to present compact analytic formulas for these 
trapped and passing guiding-center orbits expressed in 
terms of the Jacobi elliptic functions [ll|, [l2j ■ These 
compact formulas are readily amenable to explicit calcu- 
lations that do not require using simplifying limits such 
as the deeply-trapped limit (in which case the Jacobi el- 
liptic functions become simple trigonometric functions). 

The remainder of the paper is organized as follows. In 
Sec. [TT1 we introduce the simple axisymmetric tokamak 
geometry used here and in previous works @. In Sec. IIIIl 
we introduce the bounce and transit actions Ji = (Jb, Jt) 
associated with the trapped-particle (£ — b) and passing- 
particle {I — t) guiding-center orbits and their associated 
orbital frequencies uj£ = (dJt/d£)~ l , which are given 
in terms of complete elliptic integrals @. In Sec. IIV1 
we present compact expressions for the poloidal angle d 
in terms of Jacobi elliptic functions depending on the 



bounce angle Cb (for trapped-particle orbits) and the 
transit angle Ct (for passing-particle orbits). In Sec. fVl 
we introduce the canonical parallel coordinates s(Q,Ji) 
and p\\(C,i,Ji) for trapped-particle and passing-particle 
guiding-center orbits, which satisfy a canonical condition 
{s, = 1. This canonical conditions forms the basis 
of the bounce-center phase-space transformation [l3j in 
general magnetic geometry. The fact that the parallel 
coordinates explicitly satisfy the canonical conditions is 
a major result of the present work. Lastly, in Sec. lVIl the 
slow toroidal-drift guiding-center dynamics is described 
in terms of compact expressions involving the Jacobi el- 
liptic functions and our work is summarized in Sec. I VIII 



II. SIMPLE TOKAMAK GEOMETRY 

In circular, large-aspect-ratio tokamak geometry [9[, 
the magnetic-field strength is approximated as 



BftM) ~ B (l-^ cos??) 



= B e 



2B a e skT 



where the small inverse aspect ratio 

e{i>) = r(i>)/R < 1 



(1) 



(2) 



is defined in terms of the minor radius r(ip) of a circular 
magnetic surface (labeled by tp), and the poloidal angle d 
is measured from the outside equatorial plane. In Eq. ([1} , 
Bq denotes the magnetic-field strength on the magnetic 
axis (at r = 0) and B e = Bo (1 — e) denotes the magnetic- 
field strength on the (outside) equatorial plane. Lastly, 
we note that the magnitude of the toroidal magnetic is 
Btoi — B while the magnitude of the poloidal magnetic 
field at the minor radius r is B po \ ~ e(ip) B Q /q(ip), where 
q(-0) is the safety factor Q. In all expressions presented 
here, terms of order e 2 are systematically omitted (the 
symbol ~ is used whenever this approximation is first 
introduced). 

On a single magnetic surface (at fixed e), the magni- 
tude of the parallel momentum \p« |($) = ^/2m (£ — fj,B) 
of a guiding-center of mass m (at fixed total energy £, 
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magnetic moment fx, and magnetic flux tp) is 



III. BOUNCE AND TRANSIT ACTIONS 



2e uB \ . o 

sin — 

£- flB e J 2 



P||c \/l - sin 2 - = mi?,, |t?|. (3) 



Here, the equatorial parallel momentum (on the low-field 
side at i? = 0) of the guiding-center orbit is defined as 



At the lowest order on the bounce-transit time scale, 
the magnetic flux ip is frozen and thus the minor ra- 
dius r = eR and the connection length R» are con- 
stant. Hence, the parallel momentum (J3J depends only 
on the poloidal angle -d as a guiding-center moves along 
a magnetic-field line s(i?) (projected onto the poloidal 
plane at fixed toroidal angle <j>). 



P\\ c = \/2m (£ — fj,B e ) = 2y/iz mR\\u}\\, (4) 
where the bounce-transit parameter [3| is defined as 



A. Bounce action and bounce frequency for 
trapped particles 



and £ c = p\\ c /p denotes the equatorial pitch-angle co- 
ordinate. In Eq. we also defined the characteristic 
parallel frequency 



J_ / » B ° 
R\\\ e m ' 



(6) 



and the connection length R» = q(ijj) R. Using these 
definitions and Eq. the equatorial poloidal angular 
velocity is defined as 



$ c = 2i/kwi. 



(7) 



Note that the parallel frequency (0 is small compared to 
the gyrofrequency £Iq (evaluated at the magnetic axis): 



By using the infinitesimal parallel-length element ds ~ 
i?H d-d, we easily obtain the following expression for the 
bounce action for trapped particles (re < 1): 
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J b = — f p\\ ds 



Pile R\\ 



1 - k- 1 sin 2 . (11) 

2 7T 



After making the substitution sin$/2 = ^/resmyj in 
Eq. ((TTJ>, we obtain @ 



D cos 2 ip dip 

Jb = p\\ c R\\ / , 

n Jo v 1 — k si 



E(re) 



sin 2 </? 



1 - re K 



,(12) 



1 



« 1, 
2 En ' 



(8) 



where po = \/2/iBo/( m ^o) denotes the gyroradius of a 
guiding-center (with magnetic moment /x) evaluated at 
the magnetic axis. 

The condition re < 1 in Eq. Q implies that a trapped 
particle bounces back and forth between the bounce 
poloidal angles ± i?b , where 



where we wrote cos 2 ip = kT 1 [(1 — re sin 2 if) — (1 — re)] 
while K(re) and E(re) denote the complete elliptic integrals 
of the first and second kind, respectively [Til, [l^. 

By using standard properties of these complete elliptic 
integrals |11| , the bounce frequency for trapped particles 
Wb = {dJb/d£)~ 1 is expressed as 



2K(k) ; 



(13) 



2 arcsin \/re. 



(9) 



while a passing particle is in transit for k > 1 in 
Eq. ©, with the minimum parallel momentum pi| m i n = 

P|i e \/l — < Pile attained at i? = 7T. Using the defi- 
nition ([5]), the bounce-transit boundary re = 1 yields the 
trapping condition 

C 2 < 2e/(l + e) ~ 2e (10) 

in the simple tokamak magnetic field (JT|). 



where we used dn/d£ = k/(£ — /iB e ) and the relation 

El 



d_ 

dre 



E(re) - (1 - re) K(re) 



= 5 K ( K )' ( 14 ) 



In the deeply-trapped approximation, where the equato- 
rial pitch-angle coordinate is £ 2 <C 2e (so that re <C 1), 
we use K(re) ~ ir/2 and E(re) — (1 — re) K(re) ~ 7rre/4, 
and Eqs. (Tn^ - pU)) become Jb — £ £ 2 /w||o an d w b — w ||o, 
where cj|| = R7 1 £ /m. 
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FIG. 1: Normalized action J(«) = J /(mR^oj\\), given by 
Eq. (|17[) . versus the bounce-transit parameter n. The normal- 
ized bounce action (|20|1 is shown for k < 1 while the transit 
action (|21|1 is shown for k > 1. 



B. Transit action and transit frequency for passing 
particles 

For a passing particle (k > 1), the transit action is 
J t = p\\ c R\\ 



1 - k^ 1 sin 2 



2 27 



m R» to 



I w ll 



4^ 



E(0 



(15) 



The corresponding transit frequency uj t = (dJt/d£) 1 is 
expressed as 



(16) 



K(k-!) ~ K(k) ' 
where we used the relation [TTj] 

In the large-K (energetic-passing) limit, we recover J t ~ 
P|| e i?|| and wt — p|| c /(mi?||) from Eqs. (fT5]) - (fT6|) . which 
both scale as \J~rL for n^> 1. 

C. Unified representation for action integrals 

We note that the bounce action (|12l) and the transit 
action (1151) can both be combined into a single expres- 
sion for the normalized action J(k) = J/(mR 2 u)\\) for 
arbitrary bounce-transit parameter k: 



J(k) = — a{n) k 



K(«) 



cn 2 (u\n)du, (17) 



where a(n) = 1 (k < 1) or q;(k) = 1/2 (n > 1). Note that 
the factor a = 1/2 yields a consistent relation between 



FIG. 2: Normalized frequency U = w/w|| versus the bounce- 
transit parameter k. The bounce frequency (|13|) is shown for 
k < 1 while the transit frequency (|16|) is shown for ft > 1. 



the poloidal angle # and the transit angle Ct for passing 
particles [see Eq. (|2"5]) ]. 

For k < 1, we use the definition [llj for the elliptic 
integral 



,K(k) 

E(/c) = / dn (u | k) dit 
j o 



(18) 



and the identities 
cii 2 (u\k) + sn 2 (u|«:) = 1 = dn 2 (w|K) + KSn 2 (it|«) 

to obtain 

- 8 [ k{r) 2 

Jb(re) = — \k — 1 + dn (it | «)] <iw 

n Jo 1 



(19) 



E(k) - (I-k)K(k) 



(20) 



so that Jb{n = 8/7T. For k > 1, we use the conver- 
sion relation [THEH cn(u|«) = d^T/nulK -1 ) associated 
with k > 1 — > k^ 1 < 1. By using the change of variable 
v = ^ku, Eq. (TT7l) becomes 



8^ Z-^^ 1 ) 

2tt 
4^ 



dn 2 (v | k cfo 



(21) 



so that J t (K = 1) = 4/-7T. 

Figures Q] and [5] show the normalized action J = 
J I (m_R 2 oj||), defined by Eq. (|17l) . and the normalized fre- 
quency 57 = w/wii versus the bounce-transit parameter k. 
Figure [T] shows the discontinuity of the bounce and tran- 
sit actions at n = 1. This discontinuity is related to the 
fact the transit period for a marginally-passing particle 
is associated with the motion from $ = — ir to + tt, while 
the bounce period for a marginally-trapped particle is as- 
sociated with the motion from i9 = — ir to + ir and then 
back to — — ir (i.e., doubling the path taken to com- 
plete a bounce period). The bounce frequency (|13|) and 
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the transit frequency ([TBI both vanish (i.e., the corre- 
sponding periods become infinite) at the bounce-transit 
boundary k = 1 since K(ie) — > oo as x — > 1. 



IV. BOUNCE AND TRANSIT ANGLES 



The angles canonically conjugate to the bounce ac- 
tion Jt, and the transit action J t are the bounce angle 
< Cb < 27r and the transit angle — it < Ct < re- 
spectively. These angles have been expressed previously 
in terms of the incomplete elliptic integrals, which of- 
ten prevented immediate analytical applications. In the 
present Section, we use the Jacobi elliptic functions to 
give compact expressions of the poloidal angle -d in terms 
of the bounce or transit angles. 



A. Bounce angle for trapped particles 



FIG. 3: Normalized poloidal angle fl/ir versus normalized 
bounce angle Cb/ 71 " f° r ft = 0.1 (solid curve), 0.8 (dashed 
curve), and 0.999 (dotted curve). In the deeply-trapped limit 
(k <C 1), we find $ ~ $b cos Cb and $ = at Cb = t/2 for all 
values of n < 1. 



The bounce angle Cb for trapped particles (n < 1) is 
defined as 



■9b 
it/2 



3tt 

T 



7T 

2K0O J-S y/l 

7T 



1 - K-isin 2 i?'/2 
dip 



2K(k) 



re sin (/? 
mi 1 (sine | k) , (22) 



where sin 8 = sin(i9/2)/i/re and we defined Cb = tt at 
?9 = -i? b (i.e., at 8 = -tt/2). Equation (|^J can be 
inverted to give 



. 
sm — = 
2 



/^cd(2K(K)Cb/7T 



(23) 



where we used the standard notation pq = p n/qn (with p 
and q either c, s, or d) and the identities [1JJ] sn(w— 3 K) = 
sn(u + K) = cd(u), with u = 2K(/t) Cb/i" and we used the 
4 K-periodicity of sn. 

Equation (|2"31 yields the poloidal angular velocity 



(24) 



i? = - ?9 C v 7 ! - ^ sd ^2K(ft) CbA «) 



. 2^ 

sm -, 



where $ is defined in Eq. and we used the relations 
(H2]). Figure [3] shows the normalized poloidal angle fl/ir 
versus the normalized bounce angle Cb/^r for various val- 
ues of k < 1. In the standard deeply-trapped approxi- 
mation •: k 1 i (i i) ~ i)., cosCb is shown by the solid 



B. Transit angle for passing particles 

Next, the transit angle Ct for passing particles (k > 1) 
is defined as 



2 J* U W y/i-k" 1 sin 2 i?'/2 

tt r^ 2 dtp 

K(k x ) Jo yl — sin 2 93 

7? 



— — sn sm - 

K(k _1 ) V 2 



which can be easily inverted to give 



sm — = sn K(k J — 

2 \ 7T 



(25) 



/ 7c sn ( K(k) 



Ct 



(26) 



with ■& = 7r when Ct = tt (for all values of k > 1) 
and we used sn( v / Ku|K _1 ) = y/H sn(it|re) and K(k _1 ) = 
v / kK(k). We note that the definition (f25|) yields the cor- 
rect poloidal angular velocity 



= -d e cn 



(K(«Kt/7T 



= tie 



1 - K~ 



sm 




2' 

which is consistent with the definition of the transit ac- 
tion (fTSj) . Figure 0] shows the normalized poloidal angle 
■d/n versus normalized transit angle Ct for various val- 
ues of k^ 1 < 1. Note that in the energetic-passing limit 
<C 1, we find [6j ~ Ct (shown by a solid line in 
Fig. n. 



5 




0.2 0.4 0.6 0.8 1.0 



FIG. 4: Normalized poloidal angle #/7T versus normalized 
transit angle for k — 5 (solid curve), 1.05 (dashed curve), 
and 1.0001 (dotted curve). In the energetic-passing limit 
(k -1 <§C 1), we find •& ~ £t (solid curve). 
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FIG. 5: Normalized parallel momentum py /(m.-R||W|| ) versus 
the normalized bounce angle (b/^ for « = 0.1 (solid curve), 
0.8 (dashed curve), and 0.999 (dotted curve). In the deeply- 
trapped limit (k <C 1), py — —pile sm Cb (solid curve). 



V. CANONICAL PARALLEL COORDINATES 

In this Section, we derive suitable expressions for the 
canonical parallel coordinates (s,p\\), which are required 
to satisfy (with ds ~ Ri\ dd) the canonical condition 



and 



, ( _ ds dp\\ ds dp\\ 

l ' s " TkUJ o7i!k 

v <9% d£ 



P\\< 



d£ dx 



(28) 
= 1 



for both trapped and transit particles, where £ = 
27r% [«(k)/K(k)] for trapped-particle (a = 1) and 
passing-particle (a — 1/2) orbits. In Eq. (|28p . we also 
made use of the identity R\\0J\\ = {£ — \i B c ) / '{p\\ B yJTi), 
which follows from the definitions (J4j)- ([6j) . 



A. Parallel coordinates for trapped particles 



For a trapped particle, using Eq. 
mentum ([3]) is defined as 



, the parallel mo- 



'1 - k- 1 sin 2 tf/2 
P\\ c VT~k sd (xb k) • (29) 



Figure [5] shows the normalized parallel momentum 
P||/(mi?||a;||) versus the normalized bounce angle £b/"" 
for various values of n < 1, where the standard deeply- 
trapped approximation pii (£b) ~ — p\\ e sin£b is shown as 
a solid curve. 

Next, we verify that the parallel coordinates i)(£,Xb) 
and P\\(£, Xb), defined by Eqs. (|2"3")1 and (|2"9")l . satisfy the 
canonical condition (J2HJ). First, using Eq. (|23p . we obtain 



di9 
dXb 



2 v / k(1 - k) sd, 



(30) 



8tf _ y/ K /(l-K) 

~d£ ~ {£-fiB c ) 



cn + dn 2k 



<9cd 
8k 



(31) 



where the expression for ded/dn is not needed in what 
follows [II]- Next, using Eq. ([2"5]l. we obtain 



dp\\ 
dXb 



and 



d£i = -Ph 
d£ 2{£-nB e ) 



- piu y/l — k cd nd, 

(1 - 2k) , „ dsd 

v . ; sd + 2 K VI - K — 

VI — K <vAb 



(32) 



(33) 



By combining Eqs. ([30 ]) - ([33} into Eq. (|2"8]l. we obtain 
{», P||}b = ^ [« ((l-«)sd 2 + cd 2 ) ] = 1, (34) 

where we used the identity (1 — n) sd 2 + cd 2 = 1 obtained 
from Eq. (fT5|) . 



B. Parallel coordinates for passing particles 

For a passing particle, using Eq. (j2"o) . the parallel mo- 
mentum ([3]) is 



= p|| c cn (k(k) ^ ac^ , 



(35) 



where p|| m i n = p\\ e y/l — k 1 when £ t = i"- The modula- 
tion in parallel momentum 



Ap,| = P|| C (l - 

= 2mi?||W|| fy^ - v'it 7 !) (36) 
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FIG. 6: Normalized parallel momentum p||/(mi2||W||) versus 
the normalized transit angle Ct/ 71 " f° r K = 1-5 (solid curve), 
1.05 (dashed curve), and 1.0001 (dotted curve). The modula- 
tion of the parallel momentum of a passing particle decreases 



as k increases. 



decreases with increasing values of k. Figure [5] shows the 
normalized parallel momentum pu/(mR\\U}u) versus the 
normalized transit angle Ct for various values of k > 1 . 

We now verify that the parallel coordinates §{£,xt) 
and m(£ , Xt), defined by Eqs. (j2"rjj) and (pj)) . satisfy the 
canonical condition (|28p. First, using Eq. (|26p. we obtain 



B6 
dxt 



and 



d£ 



{S-iiB, 
Next, using Eq. (|35j) . we obtain 



2 yfn . 



sd + nd | 2 k — - 

OK 



d P\\ 
dxt 



- Pl|e d n sn, 



(37) 



(38) 



(39) 



and 



dp\\ 
d£ 



P\\e 



2(8 -pB e ) 
By combining Eqs. ([37 ]) -([l0" l) into Eq. 



2 k 



den 



Ok 

, we obtain 



(40) 



d_ 

dn 



{s, p\\h = 

which follows from Eq. ([T 



= 1, 



(41) 




FIG. 7: Phase portrait in (#, #)-plane of guiding-center 
trapped and transit orbits. The separatrix (dashed curves) 
separates trapped-particle orbits (inside) from the passing- 
particle orbits (outside). 



Figure [7] shows the phase portrait of the guiding-center 
trapped-particle and passing-particle orbits in the #)- 
plane, which combines expressions for the the poloidal 
angle ?9 and the normalized poloidal angular velocity 
■i?/i? c for trapped and passing particles. Here, we eas- 
ily recognize the well-known separatrix structure of the 
standard pendulum problem, which separates trapped- 
particle orbits (inside) from passing-particle orbits (out- 
side). In Sec. lVIl we will project the trapped-particle and 
passing-particle guiding-center orbits onto the poloidal 
plane (at a constant toroidal angle <$). 



VI. TOROIDAL GUIDING-CENTER 
DYNAMICS 

In this Section, we discuss the guiding-center dynam- 
ics in the toroidal direction (i.e., the direction of axisym- 
metry) and its implications on the canonically-conjugate 
drift action. So far, the toroidal angle (j> was considered 
fixed, which enabled a study of the trapped and passing 
guiding-center orbits projected onto the poloidal plane, 
and the magnetic flux ip was assumed to be frozen on the 
bounce-transit time scale. 

The axisymmetry of the tokamak magnetic geometry 
implies that the toroidal canonical momentum (or drift 
action) 



P0 = - ~ ^ 
c 

is conserved, where 
stant. 



n 



Bcfs 



- ^* 
c 



(42) 



Bq R is approximately con- 



C. Phase portrait 

We have demonstrated in Eqs. ([M]) and (|4"Tj) that the 
parallel coordinates (s,p\\), expressed in terms of Jacobi 
elliptic functions, are valid canonical coordinates for the 
trapped-particle and passing-particle guiding-center or- 
bits. 



A. Trapped-particle guiding-center orbits 

For a trapped-particle guiding-center orbit, this con- 
servation law is expressed as ip* = ipo in terms of the 
magnetic flux at the bounce points, so that Eq. (|42|) 
becomes 

ip - p\\ B^ = ip , (43) 



FIG. 8: Trapped (banana) guiding-center orbits in the nor- 
malized poloidal (x, z)-plane for n = (0.1,0.5,0.99), with the 
invariant ip* — ipo chosen as the magnetic flux ipo of the 
bounce points (shown as a dashed circle). 



FIG. 9: Passing guiding-center orbits in the normalized 
poloidal (x, z)-plane, with the invariant ip* = ipt chosen as 
the magnetic flux ipt on the inside equatorial plane (shown 
as a dashed circle). The dotted and solid circles are passing 
guiding-center orbits with k = 1.05 and 1.5, respectively. 



where pu = pu/(mQ,). By substituting the Taylor ap- 
proximation ip ~ ip + A r (e B^/qq), the normalized 
radial deviation Aeb = e — £o (from the magnetic surface 
ip = ipo) is expressed as 



Ae b 



eo R 



2 S y/n (1 - k) sd (2K(k) CbA k) 



(44) 

where, using Eq. (|8]l. we have defined the trapped-orbit 
parameter 



<7o Po 



rfa R 



(45) 



The maximum radial deviation (|4"4")l is (Aeb) max = 
25 2(5. 

A trapped-particle guiding-center orbit can be pro- 
jected onto the normalized poloidal plane (x = x/R,z = 
z/R) and represented as 



s(Cb,«) = 1 + (e + Ae b ) cos?9(Cb,K) 
s(Cb,«) = (eo + Ae b ) sini9(C b ,K) 



(46) 



Figure [5] shows the trapped (banana) guiding-center or- 
bits for (eo,(5) = (0.5,0.1) for various values of « < 1, 
with the magnetic surface ip = ipo shown as a dashed 
circle. We note that, in the simple orbit topology consid- 
ered here, a trapped-particle orbit does not enclose the 
magnetic axis [15fl . 



magnetic flux ip* = ip t at the inside equatorial mid- 
plane (at i? = tt). By using the Taylor approximation 
ip ~ ipt + Ar (e t B^/qt), the normalized radial devia- 
tion Ae t = e — e t (from the magnetic surface ip = ipt) is 
expressed as 



Ae t 



2(5 



v^dn K(k" 1 ) 



Ct 



(47) 

where orbit parameter (|45j) is now evaluated at 0t (note 
that Aet vanishes on the inside equatorial plane where 
P\\ = Pllmin = Pile Vl — k _1 ). The maximum radial de- 
viation (|47l) is (Ae t )max = 2(5 (-/k — v 7 ^ — 1), i-C, the 
largest radial deviation for a passing-particle orbit oc- 
curs for marginally-passing particles (« — 1 <C 1). 

A passing-particle guiding-center orbit can be pro- 
jected onto the normalized poloidal plane and repre- 
sented as 



x((t,K) = 1 + (e t + Ae t ) costf(Ct,«0 
z((t,K) = (e t + Ae t ) sinz?(Ct,K) 



(48) 



Figure [9] shows the passing guiding-center orbits for 
(e t ,<5) = (0.5,0.1) for various values of k > 1, with the 
magnetic surface ip = ipt shown as a dashed circle. 



C. Bounce-averaged toroidal drift frequency 



B. Passing-particle guiding-center orbits The bounce- averaged toroidal drift frequency [|| 

For a passing-particle guiding-center orbit, the con- / • k\ — C P / 



servation law (l4lfl) may be expressed in terms of the U>d y ^ / e \ dip ' ' 
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where we used Eq. (fl"4]) and 8k/ de = (1/2 — n)/e. The 
toroidal drift frequency (|50j) can thus be approximately 
expressed as 



0_ 
~de 



,R 2 



(52) 



Hence, the toroidal drift reversal occurs when the radial 
derivative of the bounce action (fl~2]) changes sign. 



FIG. 10: Normalized drift frequency [E(«)/K(k)] - 1/2 as a 
function of k. Note that the toroidal drift frequency (|50|) 
changes sign for k > k t ~ 0.826. 



VII. SUMMARY 



is approximated (ignoring magnetic shear) as 



2 f iB a 



mO ei? 2 VK(k) 



E(«) 



where we used 8B/d4> — ~(l/ e R 2 ) cos-d. We note 
that the bounce-averaged toroidal drift frequency (|50j) 
exhibits a drift reversal for k > n r ~ 0.826, where 
K(k t ) = 2 E(« r ) (see Fig.[TDl) ; a more complete expression 
for the bounce-averaged toroidal drift frequency can be 
found elsewhere [l|. 

Lastly, the toroidal drift frequency (|50|) can also be 
expressed in terms of a radial derivative of the bounce 
action as follows. First, we write 



de 



4m o 
^11 

7T 6 11 



1 



K 



i-,|K 



yB 



The analytical formulas presented in this paper pro- 
vide compact expressions in terms of the Jacobi ellip- 
tic functions (cn, sn, dn) that describe trapped-particle 

(50) ( K < 1) an d passing-particle (k > 1) guiding-center or- 
bits in simple axisymmetric tokamak geometry for ar- 
bitrary bounce-transit parameter k. The parallel coordi- 
nates (s,p||) derived in Sec. El which satisfy the canonical 
condition for trapped particles and (|4"Tj) for passing 
particles, guarantee that the bounce-center transforma- 
tion Q in axisymmetric tokamak geometry can be explic- 
itly carried out beyond the deeply-trapped and energetic- 
passing limits @. 

Applications of the compact Jacobi-elliptic represen- 
tations of trapped-particle and passing-particle guiding- 
center orbits include the derivation of bounce-center- 
averaged fluctuating potentials in bounce-kinetic theory 
[fllll and the derivation of a bounce-center Fokker-Planck 
collision operator @, [1[ beyond the zero-banana-width 

(51) approximation. 
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